<html>
 <head>
  <meta charset="utf-8"/>
  <meta content="width=device-width, initial-scale=1, maximum-scale=1, user-scalable=no" name="viewport"/>
  <title>
   谢益辉的回复  | 数螺 | NAUT IDEA
  </title>
  <link href="http://cdn.bootcss.com/bootstrap/3.3.6/css/bootstrap-theme.min.css" rel="stylesheet"/>
  <link href="http://cdn.bootcss.com/bootstrap/3.3.6/css/bootstrap.min.css" rel="stylesheet"/>
  <style type="text/css">
   #xmain img {
                  max-width: 100%;
                  display: block;
                  margin-top: 10px;
                  margin-bottom: 10px;
                }

                #xmain p {
                    line-height:150%;
                    font-size: 16px;
                    margin-top: 20px;
                }

                #xmain h2 {
                    font-size: 24px;
                }

                #xmain h3 {
                    font-size: 20px;
                }

                #xmain h4 {
                    font-size: 18px;
                }


                .header {
	           background-color: #0099ff;
	           color: #ffffff;
	           margin-bottom: 20px;
	        }

	        .header p {
                  margin: 0px;
                  padding: 10px 0;
                  display: inline-block;  
                  vertical-align: middle;
                  font-size: 16px;
               }

               .header a {
                 color: white;
               }

              .header img {
                 height: 25px;
              }
  </style>
  <script src="http://cdn.bootcss.com/jquery/3.0.0/jquery.min.js">
  </script>
  <script src="http://cdn.mathjax.org/mathjax/latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML" type="text/javascript">
   MathJax.Hub.Config({elements: ["bbpress-forums"]});
  </script>
  <script src="http://nautstatic-10007657.file.myqcloud.com/static/css/readability.min.js" type="text/javascript">
  </script>
  <script type="text/javascript">
   $(document).ready(function() {
                 var loc = document.location;
                 var uri = {
                  spec: "http://cos.name/cn/profile/1/replies/page/4/",
                  host: "http://cos.name",
                  prePath: "http://cos.name",
                  scheme: "http",
                  pathBase: "http://cos.name/"
                 };
    
                 var documentClone = document.cloneNode(true);
                 var article = new Readability(uri, documentClone).parse();
     
                 document.getElementById("xmain").innerHTML = article.content;
                });
  </script>
  <!-- 1466443423: Accept with keywords: (title(0.0):论坛, topn(0.4):暴力,方差,办法,帖子,分布,字体,用户名,意义,结果,分位,参考文献,可视化,随机数,时间序列,方法,正态分布,中文,函数,正则表达式,公式,回归,楼主,问题,篇文章,随机变量,有时候,主题,榆木脑袋,作者,文献).-->
 </head>
 <body class="bbp-user-page single singular bbpress paged-4 single-author sidebar" onload="">
  <div class="header">
   <div class="container">
    <div class="row">
     <div class="col-xs-6 col-sm-6 text-left">
      <a href="/databee">
       <img src="http://nautidea-10007657.cos.myqcloud.com/logo_white.png"/>
      </a>
      <a href="/databee">
       <p>
        数螺
       </p>
      </a>
     </div>
     <div class="hidden-xs col-sm-6 text-right">
      <p>
       致力于数据科学的推广和知识传播
      </p>
     </div>
    </div>
   </div>
  </div>
  <div class="container text-center">
   <h1>
    谢益辉的回复
   </h1>
  </div>
  <div class="container" id="xmain">
   <div class="hfeed site" id="page">
    <header class="site-header" id="masthead" role="banner">
     <div id="cos-logo">
      <a href="http://cos.name/cn">
       <img src="http://cos.name/cn/wp-content/themes/COS-forest/images/headers/cos-logo.png"/>
      </a>
     </div>
     <div class="navbar" id="navbar">
      <nav class="navigation main-navigation" id="site-navigation" role="navigation">
       <h3 class="menu-toggle">
        菜单
       </h3>
       <div class="menu-%e8%8f%9c%e5%8d%951-container">
        <ul class="nav-menu" id="menu-%e8%8f%9c%e5%8d%951">
         <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-home menu-item-407772" id="menu-item-407772">
          <a href="http://cos.name/cn/">
           论坛首页
          </a>
         </li>
         <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-407773" id="menu-item-407773">
          <a href="http://cos.name/cn/forums/">
           讨论区
          </a>
         </li>
         <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-407774" id="menu-item-407774">
          <a href="http://cos.name/cn/wp-login.php?action=register">
           注册
          </a>
         </li>
         <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-407819" id="menu-item-407819">
          <a href="http://cos.name/">
           主站
          </a>
         </li>
        </ul>
       </div>
      </nav>
      <!-- #site-navigation -->
     </div>
     <!-- #navbar -->
    </header>
    <!-- #masthead -->
    <div class="site-main" id="main">
     <div class="content-area" id="primary">
      <div class="site-content" id="content" role="main">
       <article class="post-0 type- status-publish hentry" id="post-0">
        <header class="entry-header">
         <h1 class="entry-title">
          谢益辉
         </h1>
        </header>
        <!-- .entry-header -->
        <div class="entry-content">
         <div id="bbpress-forums">
          <div id="bbp-user-wrapper">
           <div id="bbp-single-user-details">
            <div id="bbp-user-avatar">
             <span class="vcard">
              <a class="url fn n" href="http://cos.name/cn/profile/1/" rel="me" title="谢益辉">
               <img src="http://sdn.geekzu.org/avatar/1022d8e6ebc94e8f6bca9a86cebe312a?s=150&amp;d=monsterid&amp;r=g"/>
              </a>
             </span>
            </div>
            <!-- #author-avatar -->
            <div id="bbp-user-navigation">
             <ul>
              <li class="">
               <span class="vcard bbp-user-profile-link">
                <a class="url fn n" href="http://cos.name/cn/profile/1/" rel="me" title="谢益辉的档案">
                 档案
                </a>
               </span>
              </li>
              <li class="">
               <span class="bbp-user-topics-created-link">
                <a href="http://cos.name/cn/profile/1/topics/" title="谢益辉发起的主题">
                 发起的主题
                </a>
               </span>
              </li>
              <li class="current">
               <span class="bbp-user-replies-created-link">
                <a href="http://cos.name/cn/profile/1/replies/" title="谢益辉创建的回复">
                 创建的回复
                </a>
               </span>
              </li>
              <li class="">
               <span class="bbp-user-favorites-link">
                <a href="http://cos.name/cn/profile/1/favorites/" title="谢益辉的收藏夹">
                 收藏夹
                </a>
               </span>
              </li>
             </ul>
            </div>
            <!-- #bbp-user-navigation -->
           </div>
           <!-- #bbp-single-user-details -->
           <div id="bbp-user-body">
            <div class="bbp-user-replies-created" id="bbp-user-replies-created">
             <h2 class="entry-title">
              回复的主题
             </h2>
             <div class="bbp-user-section">
              <div class="bbp-pagination">
               <div class="bbp-pagination-count">
                查看 15 个帖子 - 46 到 60（总计 8,682 个）
               </div>
               <div class="bbp-pagination-links">
                <a class="prev page-numbers" href="http://cos.name/cn/profile/1/replies/page/3/">
                 ←
                </a>
                <a class="page-numbers" href="http://cos.name/cn/profile/1/replies/">
                 1
                </a>
                <span class="page-numbers dots">
                 …
                </span>
                <a class="page-numbers" href="http://cos.name/cn/profile/1/replies/page/3/">
                 3
                </a>
                <span class="page-numbers current">
                 4
                </span>
                <a class="page-numbers" href="http://cos.name/cn/profile/1/replies/page/5/">
                 5
                </a>
                <span class="page-numbers dots">
                 …
                </span>
                <a class="page-numbers" href="http://cos.name/cn/profile/1/replies/page/579/">
                 579
                </a>
                <a class="next page-numbers" href="http://cos.name/cn/profile/1/replies/page/5/">
                 →
                </a>
               </div>
              </div>
              <ul class="forums bbp-replies" id="topic-0-replies">
               <li class="bbp-header">
                <div class="bbp-reply-author">
                 作者
                </div>
                <!-- .bbp-reply-author -->
                <div class="bbp-reply-content">
                 帖子
                </div>
                <!-- .bbp-reply-content -->
               </li>
               <!-- .bbp-header -->
               <li class="bbp-body">
                <div class="bbp-reply-header" id="post-392056">
                 <div class="bbp-meta">
                  <span class="bbp-reply-post-date">
                   2014年1月6日 下午9:29
                  </span>
                  <span class="bbp-header">
                   回复：
                   <a class="bbp-topic-permalink" href="http://cos.name/cn/topic/147769/">
                    R能不能绘制这种图形？？
                   </a>
                  </span>
                  <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/147769/page/2/#post-392056">
                   23 楼
                  </a>
                  <span class="bbp-admin-links">
                  </span>
                 </div>
                 <!-- .bbp-meta -->
                </div>
                <!-- #post-392056 -->
                <div class="odd bbp-parent-forum-999 bbp-parent-topic-147769 bbp-reply-position-23 user-id-1 post-392056 reply type-reply status-publish hentry">
                 <div class="bbp-reply-author">
                  <a class="bbp-author-avatar" href="http://cos.name/cn/profile/1/" rel="nofollow" title="查看谢益辉的档案">
                   <img src="http://sdn.geekzu.org/avatar/1022d8e6ebc94e8f6bca9a86cebe312a?s=80&amp;d=monsterid&amp;r=g"/>
                  </a>
                  <br/>
                  <a class="bbp-author-name" href="http://cos.name/cn/profile/1/" rel="nofollow" title="查看谢益辉的档案">
                   谢益辉
                  </a>
                  <br/>
                  <div class="bbp-author-role">
                   站长
                  </div>
                 </div>
                 <!-- .bbp-reply-author -->
                 <div class="bbp-reply-content">
                  <p>
                   完全颠覆了我们对字体的印象，字体本来是字符到坐标信息的对应，而这个坐标信息是可以随意给的，没有人规定A的字体必须画出A的形状。
                  </p>
                 </div>
                 <!-- .bbp-reply-content -->
                </div>
                <!-- .reply -->
                <div class="bbp-reply-header" id="post-389896">
                 <div class="bbp-meta">
                  <span class="bbp-reply-post-date">
                   2013年12月28日 上午7:44
                  </span>
                  <span class="bbp-header">
                   回复：
                   <a class="bbp-topic-permalink" href="http://cos.name/cn/topic/142002/">
                    多重比较可视化求助
                   </a>
                  </span>
                  <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/142002/page/2/#post-389896">
                   19 楼
                  </a>
                  <span class="bbp-admin-links">
                  </span>
                 </div>
                 <!-- .bbp-meta -->
                </div>
                <!-- #post-389896 -->
                <div class="even bbp-parent-forum-999 bbp-parent-topic-142002 bbp-reply-position-19 user-id-1 post-389896 reply type-reply status-publish hentry">
                 <div class="bbp-reply-author">
                  <a class="bbp-author-avatar" href="http://cos.name/cn/profile/1/" rel="nofollow" title="查看谢益辉的档案">
                   <img src="http://sdn.geekzu.org/avatar/1022d8e6ebc94e8f6bca9a86cebe312a?s=80&amp;d=monsterid&amp;r=g"/>
                  </a>
                  <br/>
                  <a class="bbp-author-name" href="http://cos.name/cn/profile/1/" rel="nofollow" title="查看谢益辉的档案">
                   谢益辉
                  </a>
                  <br/>
                  <div class="bbp-author-role">
                   站长
                  </div>
                 </div>
                 <!-- .bbp-reply-author -->
                 <div class="bbp-reply-content">
                  <p>
                   或许的确是一种更保守的办法（SMM方法本身已经是保守的办法了），但我想不出它在什么情况下有用，主要问题是这样：一共有r组数据，要做r(r-1)次比较，这个方法的前提是每一组均值的差异底下除的是一些相互独立的方差的估计，这个条件难以满足（当然不是说不能造出来一些独立的估计，只是多组混合在一起的方差估计会更好一些）。话说回来呢，SMM方法本身已经是“犯规”了，因为那些分子|Xi – Xj|之间显然也不独立，这也是它的保守性所在。五十步已经有些离谱了，有没有必要一百步呢？我估计这里面已经没什么文章可做了。要不是你开这个话题，其实我都不知道有SMM这么个东西的存在，教科书上都是教的Tukey’s HSD（基于Studentized Range而不是Modulus）。
                  </p>
                 </div>
                 <!-- .bbp-reply-content -->
                </div>
                <!-- .reply -->
                <div class="bbp-reply-header" id="post-389611">
                 <div class="bbp-meta">
                  <span class="bbp-reply-post-date">
                   2013年12月27日 下午5:41
                  </span>
                  <span class="bbp-header">
                   回复：
                   <a class="bbp-topic-permalink" href="http://cos.name/cn/topic/142002/">
                    多重比较可视化求助
                   </a>
                  </span>
                  <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/142002/page/2/#post-389611">
                   16 楼
                  </a>
                  <span class="bbp-admin-links">
                  </span>
                 </div>
                 <!-- .bbp-meta -->
                </div>
                <!-- #post-389611 -->
                <div class="odd bbp-parent-forum-999 bbp-parent-topic-142002 bbp-reply-position-16 user-id-1 post-389611 reply type-reply status-publish hentry">
                 <div class="bbp-reply-author">
                  <a class="bbp-author-avatar" href="http://cos.name/cn/profile/1/" rel="nofollow" title="查看谢益辉的档案">
                   <img src="http://sdn.geekzu.org/avatar/1022d8e6ebc94e8f6bca9a86cebe312a?s=80&amp;d=monsterid&amp;r=g"/>
                  </a>
                  <br/>
                  <a class="bbp-author-name" href="http://cos.name/cn/profile/1/" rel="nofollow" title="查看谢益辉的档案">
                   谢益辉
                  </a>
                  <br/>
                  <div class="bbp-author-role">
                   站长
                  </div>
                 </div>
                 <!-- .bbp-reply-author -->
                 <div class="bbp-reply-content">
                  <p>
                   最后一次更新：三个字，我错了。
                  </p>
                  <p>
                   原来的公式是对的，你在6楼的第一种模拟比较接近正确的答案。我之所以错了是因为这个问题并不是独立的t随机变量，尽管每个随机变量的确是t分布，但不独立，因为每个正态都除以同一个卡方的平方根，分母相同，分子独立，这些商之间并不独立。昨晚躺床上想起来这么一茬事，睡得比平时都快。
                  </p>
                  <p>
                   真正的模拟应该是这样的：
                  </p>
                  <p>
                  </p>
                  <pre class="highlight ">rsmm = function(n, r, df) {
  x_max = apply(matrix(abs(rnorm(r * n)), r, n), 2, max)
  s_chi = sqrt(rchisq(n, df)/df)
  x_max/s_chi
}
qsmm = function(q, r, df, n = 100000) quantile(rsmm(n, r, df), q)</pre>
                  <p>
                  </p>
                  <p>
                   7楼的积分仍然是目前我能想到的正确的答案。这回是真丢人了，榆木脑袋宣布结案，请不要轻信所谓的科班或权威。[s:18]
                  </p>
                  <p>
                   另外一种基于模拟的方法是：
                  </p>
                  <p>
                  </p>
                  <pre class="highlight ">psmm = function(x, r, df, n = 1000000) {
  chi = sqrt(rchisq(n, df)/df)
  mean((2 * pnorm(x * chi) - 1)^r)
}
qsmm = function(q, r, df) {
  if (!is.finite(df)) return(qnorm(1 - .5 *(1 - q^(1/r))))
  res = uniroot(function(x, r, df, q) {
    psmm(x, r = r, df = df) - q
  }, c(0, 100), r = r, df = df, q = q)
  res$root
}</pre>
                  <p>
                  </p>
                  <p>
                   它主要是大数定律的应用，注意到分布函数2.2是一个积分的形式，实际上就是一个期望E((2Phi(cX) -1)^r)，其中X服从变形的Chi分布。期望可以通过随机数的实现并求函数值的平均来近似，即E(g(X))约等于mean(g(x1), …, g(xn))，x_i来自于它自己的分布。模拟的方法有不确定性，在这个问题里不如暴力积分方法。
                  </p>
                 </div>
                 <!-- .bbp-reply-content -->
                </div>
                <!-- .reply -->
                <div class="bbp-reply-header" id="post-389514">
                 <div class="bbp-meta">
                  <span class="bbp-reply-post-date">
                   2013年12月27日 上午4:16
                  </span>
                  <span class="bbp-header">
                   回复：
                   <a class="bbp-topic-permalink" href="http://cos.name/cn/topic/142002/">
                    多重比较可视化求助
                   </a>
                  </span>
                  <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/142002/#post-389514">
                   14 楼
                  </a>
                  <span class="bbp-admin-links">
                  </span>
                 </div>
                 <!-- .bbp-meta -->
                </div>
                <!-- #post-389514 -->
                <div class="even bbp-parent-forum-999 bbp-parent-topic-142002 bbp-reply-position-14 user-id-1 post-389514 reply type-reply status-publish hentry">
                 <div class="bbp-reply-author">
                  <a class="bbp-author-avatar" href="http://cos.name/cn/profile/1/" rel="nofollow" title="查看谢益辉的档案">
                   <img src="http://sdn.geekzu.org/avatar/1022d8e6ebc94e8f6bca9a86cebe312a?s=80&amp;d=monsterid&amp;r=g"/>
                  </a>
                  <br/>
                  <a class="bbp-author-name" href="http://cos.name/cn/profile/1/" rel="nofollow" title="查看谢益辉的档案">
                   谢益辉
                  </a>
                  <br/>
                  <div class="bbp-author-role">
                   站长
                  </div>
                 </div>
                 <!-- .bbp-reply-author -->
                 <div class="bbp-reply-content">
                  <p>
                   我提到的两篇文献在这里：
                   <br/>
                   <a href="https://www.dropbox.com/s/1g9xjzv87t7otdh/1268584.pdf" rel="nofollow">
                    https://www.dropbox.com/s/1g9xjzv87t7otdh/1268584.pdf
                   </a>
                   <br/>
                   <a href="https://www.dropbox.com/s/fov1m4mm3tlqv4a/2334520.pdf" rel="nofollow">
                    https://www.dropbox.com/s/fov1m4mm3tlqv4a/2334520.pdf
                   </a>
                  </p>
                  <p>
                   那本书确实是Simultaneous Statistical Inference。亚马逊上有预览，不过预览里关键的几页看不到，我们图书馆有这本书，我懒得爬出去找了。
                  </p>
                  <p>
                   今天我再想了一下这个问题，意识到它其实比我原先想象的要简单得多，因为我比较懒，能模拟的情况下先考虑模拟，但模拟其实是最劣等的方案，最理想的方案当然是找到解析解，对这个问题来说解析解不太可能。折衷一下，半计算半解析：
                  </p>
                  <p>
                   若X服从t分布，累积分布函数是F(x)，密度函数是f(x)，令Y = |X|，那么密度g(y) = 2f(y)，y &gt; 0；分布G(y) = 2F(y) – 1。对r个独立的t分布随机变量X1, …, Xn，令Z = max(|X1|, …, |Xn|)，那么Z分布H(z) = (2F(z) – 1)^r。这都是些比较初等的概率推导。
                  </p>
                  <p>
                   既然分布函数能以很简单的形式写出来，代码也就大大简化了，可以直接用t分布的分布函数
                   <code>
                    pt()
                   </code>
                   ：
                  </p>
                  <p>
                  </p>
                  <pre class="highlight ">psmm = function(x, r, nu) {
  ifelse(x &gt; 0, (2 * pt(x, nu) - 1)^r, 0)
}
# 分位数函数仍然用H(x) - q = 0求解x
qsmm = function(q, r, nu) {
  res = uniroot(function(x, r, nu, q) {
    psmm(x, r = r, nu = nu) - q
  }, c(0, 100), r = r, nu = nu, q = q)
  res$root
}
# 算0.99分位数
qsmm(.99, r=3, nu=5)
ks = 3:20
vs = c(5, 7, 10, 12, 16, 20, 24, 30, 40, 60, 120, Inf)
q0.01 = matrix(nrow = length(ks), ncol = length(vs), dimnames = list(ks, vs))
for (k in ks) {
  for (v in vs) {
    q0.01[as.character(k), as.character(v)] = qsmm(.99, r = k*(k-1)/2, nu = v)
  }
}
cbind('r*' = choose(ks, 2), round(q0.01, 3))
    r*      5     7    10    12    16    20    24    30    40    60   120   Inf
3    3  5.243 4.353 3.825 3.647 3.443 3.329 3.257 3.188 3.121 3.056 2.994 2.934
4    6  6.133 4.941 4.256 4.029 3.771 3.629 3.539 3.453 3.370 3.291 3.215 3.143
5   10  6.862 5.404 4.584 4.315 4.013 3.848 3.744 3.644 3.549 3.459 3.372 3.289
6   15  7.491 5.791 4.852 4.547 4.206 4.021 3.905 3.794 3.689 3.589 3.493 3.402
7   21  8.051 6.126 5.079 4.742 4.368 4.165 4.038 3.918 3.803 3.695 3.591 3.493
8   28  8.558 6.424 5.278 4.912 4.506 4.288 4.152 4.023 3.900 3.784 3.673 3.569
9   36  9.024 6.693 5.454 5.061 4.628 4.396 4.251 4.114 3.984 3.861 3.744 3.634
10  45  9.457 6.939 5.614 5.196 4.737 4.491 4.339 4.194 4.058 3.929 3.807 3.691
11  55  9.862 7.166 5.760 5.319 4.835 4.577 4.417 4.267 4.124 3.989 3.862 3.742
12  66 10.244 7.378 5.894 5.431 4.925 4.656 4.489 4.332 4.184 4.044 3.912 3.787
13  78 10.606 7.576 6.019 5.535 5.008 4.728 4.555 4.392 4.238 4.094 3.957 3.829
14  91 10.950 7.762 6.135 5.632 5.084 4.794 4.615 4.447 4.288 4.139 3.999 3.866
15 105 11.279 7.939 6.245 5.722 5.156 4.856 4.672 4.498 4.335 4.181 4.037 3.901
16 120 11.595 8.107 6.348 5.807 5.223 4.914 4.724 4.546 4.378 4.221 4.073 3.933
17 136 11.898 8.267 6.446 5.888 5.286 4.968 4.773 4.590 4.418 4.257 4.106 3.963
18 153 12.190 8.419 6.538 5.964 5.345 5.020 4.820 4.632 4.456 4.291 4.137 3.991
19 171 12.472 8.566 6.627 6.037 5.402 5.068 4.863 4.672 4.492 4.324 4.166 4.018
20 190 12.745 8.707 6.711 6.106 5.455 5.114 4.905 4.709 4.526 4.354 4.193 4.042</pre>
                  <p>
                  </p>
                  <p>
                   文献中的分位数表是下面这个，可以看出自由度无穷大的时候能对得上（正态分布），但自由度比较小的时候他的分位数值总体偏小（相应置信区间的长度也会小）：
                  </p>
                  <p>
                   <img src="http://i.imgur.com/uxzspO6.png"/>
                  </p>
                  <p>
                   要么是 a) 参考文献13（Stoline等人）的Fortran程序写的有问题，要么是 b) R的t分布函数写的有问题，要么是 c) 那1966年的书推导有问题。现在看来 a) 和 c) 应该是一致的，因为他们的结果我用暴力积分的办法可以重复出来。那么要么R错了，要么公式错了。希望不是R搞错了，不然感觉真的不会再爱了。
                  </p>
                  <p>
                   20分钟后更新
                   <br/>
                   ===
                  </p>
                  <p>
                   再看了一遍公式2.2，发现它真的错了耶。
                  </p>
                  <p>
                   <img src="http://i.imgur.com/JNE3Wcz.png"/>
                  </p>
                  <p>
                   那个幂r应该提到积分的外面去。我了个去，我这榆木脑袋居然能发现一处错误的公式，生平第一次哎。叩谢楼主的这份圣诞礼物。
                  </p>
                 </div>
                 <!-- .bbp-reply-content -->
                </div>
                <!-- .reply -->
                <div class="bbp-reply-header" id="post-388677">
                 <div class="bbp-meta">
                  <span class="bbp-reply-post-date">
                   2013年12月25日 上午7:37
                  </span>
                  <span class="bbp-header">
                   回复：
                   <a class="bbp-topic-permalink" href="http://cos.name/cn/topic/142002/">
                    多重比较可视化求助
                   </a>
                  </span>
                  <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/142002/#post-388677">
                   7 楼
                  </a>
                  <span class="bbp-admin-links">
                  </span>
                 </div>
                 <!-- .bbp-meta -->
                </div>
                <!-- #post-388677 -->
                <div class="odd bbp-parent-forum-999 bbp-parent-topic-142002 bbp-reply-position-7 user-id-1 post-388677 reply type-reply status-publish hentry">
                 <div class="bbp-reply-author">
                  <a class="bbp-author-avatar" href="http://cos.name/cn/profile/1/" rel="nofollow" title="查看谢益辉的档案">
                   <img src="http://sdn.geekzu.org/avatar/1022d8e6ebc94e8f6bca9a86cebe312a?s=80&amp;d=monsterid&amp;r=g"/>
                  </a>
                  <br/>
                  <a class="bbp-author-name" href="http://cos.name/cn/profile/1/" rel="nofollow" title="查看谢益辉的档案">
                   谢益辉
                  </a>
                  <br/>
                  <div class="bbp-author-role">
                   站长
                  </div>
                 </div>
                 <!-- .bbp-reply-author -->
                 <div class="bbp-reply-content">
                  <p>
                   你后面更新的是对的，我刚开始没有仔细看公式。我重新看了一下相关的文献，这个SMM确实应该是t分布的绝对值，随机数可以这样生成：
                  </p>
                  <p>
                  </p>
                  <pre class="highlight ">rsmm = function(n, m, df) apply(matrix(abs(rt(m * n, df)), m, n), 2, max)</pre>
                  <p>
                  </p>
                  <p>
                   不小心被你带到这条沟里来了，索性看了一下这篇文章：http://www.jstor.org/discover/10.2307/1268584 用他给的公式(2.2)暴力算了一下他的表1，结果能吻合：
                  </p>
                  <p>
                  </p>
                  <pre class="highlight ">psmm_x = function(x, c, r, nu) {
  snu = sqrt(nu)
  sx  = snu * x  # for the scaled Chi distribution
  lgx = log(snu) - lgamma(nu/2) + (1 - nu/2) * log(2) + (nu - 1) * log(sx) + (-sx^2/2)
  exp(r * log(2 * pnorm(c * x) - 1) + lgx)
}
psmm = function(x, r, nu) {
  res = integrate(psmm_x, 0, Inf, c = x, r = r, nu = nu)
  res$value
}
qsmm = function(q, r, nu) {
  r = (r * (r - 1)/2)
  if (!is.finite(nu)) return(qnorm(1 - .5 *(1 - q^(1/r))))
  res = uniroot(function(c, r, nu, q) {
    psmm(c, r = r, nu = nu) - q
  }, c(0, 100), r = r, nu = nu, q = q)
  res$root
}
qsmm(.99, r=3, nu=5)
qsmm(.99, r=3, nu=Inf)
ks = 3:20
vs = c(5, 7, 10, 12, 16, 20, 24, 30, 40, 60, 120, Inf)
q0.01 = matrix(nrow = length(ks), ncol = length(vs), dimnames = list(ks, vs))
for (k in ks) {
  for (v in vs) {
    q0.01[as.character(k), as.character(v)] = qsmm(.99, r = k, nu = v)
  }
}
round(q0.01, 3)</pre>
                  <p>
                  </p>
                  <p>
                   代码写得比较混乱，大概意思是用uniroot()求F(x) = q的根（即分位数），累积概率分布函数是用积分integrate()求的，被积函数为了数值稳定性先取log()再exp()。当自由度无穷大时，t分布就是正态分布，文中的公式(3.5)给出了分位数算法，就不用算积分那么麻烦了。
                  </p>
                  <p>
                   虽然结果能吻合，但还是疑云重重，因为毕竟模拟出来的结果和他给的结果还是有差异的，有时候大，有时候小。说实话我有点怀疑那个公式有问题，不过作者说是从参考文献5和8里来的，而我找到了5，发现里面没有他说的公式2.2（童话里都是骗银滴），文献8是一本1966年的书，书这种东西比较难找，我挖坑只能挖到这里了，接下来把锹还给你接着挖吧，搞不好能发篇文章正本清源发现大家用了四十多年的多重比较是错的呢。
                  </p>
                 </div>
                 <!-- .bbp-reply-content -->
                </div>
                <!-- .reply -->
                <div class="bbp-reply-header" id="post-388163">
                 <div class="bbp-meta">
                  <span class="bbp-reply-post-date">
                   2013年12月23日 上午3:52
                  </span>
                  <span class="bbp-header">
                   回复：
                   <a class="bbp-topic-permalink" href="http://cos.name/cn/topic/140990/">
                    googlevis 支持中文吗？
                   </a>
                  </span>
                  <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/140990/#post-388163">
                   2 楼
                  </a>
                  <span class="bbp-admin-links">
                  </span>
                 </div>
                 <!-- .bbp-meta -->
                </div>
                <!-- #post-388163 -->
                <div class="even bbp-parent-forum-1011 bbp-parent-topic-140990 bbp-reply-position-2 user-id-1 post-388163 reply type-reply status-publish hentry">
                 <div class="bbp-reply-author">
                  <a class="bbp-author-avatar" href="http://cos.name/cn/profile/1/" rel="nofollow" title="查看谢益辉的档案">
                   <img src="http://sdn.geekzu.org/avatar/1022d8e6ebc94e8f6bca9a86cebe312a?s=80&amp;d=monsterid&amp;r=g"/>
                  </a>
                  <br/>
                  <a class="bbp-author-name" href="http://cos.name/cn/profile/1/" rel="nofollow" title="查看谢益辉的档案">
                   谢益辉
                  </a>
                  <br/>
                  <div class="bbp-author-role">
                   站长
                  </div>
                 </div>
                 <!-- .bbp-reply-author -->
                 <div class="bbp-reply-content">
                  <p>
                   Flash中的中文从来就是让人头疼的问题，基于JS的可视化工具应该不存在这样的问题，支不支持就看开发者有没有兴趣了。
                  </p>
                 </div>
                 <!-- .bbp-reply-content -->
                </div>
                <!-- .reply -->
                <div class="bbp-reply-header" id="post-388131">
                 <div class="bbp-meta">
                  <span class="bbp-reply-post-date">
                   2013年12月22日 上午6:08
                  </span>
                  <span class="bbp-header">
                   回复：
                   <a class="bbp-topic-permalink" href="http://cos.name/cn/topic/142002/">
                    多重比较可视化求助
                   </a>
                  </span>
                  <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/142002/#post-388131">
                   5 楼
                  </a>
                  <span class="bbp-admin-links">
                  </span>
                 </div>
                 <!-- .bbp-meta -->
                </div>
                <!-- #post-388131 -->
                <div class="odd bbp-parent-forum-999 bbp-parent-topic-142002 bbp-reply-position-5 user-id-1 post-388131 reply type-reply status-publish hentry">
                 <div class="bbp-reply-author">
                  <a class="bbp-author-avatar" href="http://cos.name/cn/profile/1/" rel="nofollow" title="查看谢益辉的档案">
                   <img src="http://sdn.geekzu.org/avatar/1022d8e6ebc94e8f6bca9a86cebe312a?s=80&amp;d=monsterid&amp;r=g"/>
                  </a>
                  <br/>
                  <a class="bbp-author-name" href="http://cos.name/cn/profile/1/" rel="nofollow" title="查看谢益辉的档案">
                   谢益辉
                  </a>
                  <br/>
                  <div class="bbp-author-role">
                   站长
                  </div>
                 </div>
                 <!-- .bbp-reply-author -->
                 <div class="bbp-reply-content">
                  <p>
                   既然没有解析解，学术以及帝都雾都大，在家偷懒敲键盘整个模拟解应该可以吧。反正计算便宜，就多算点好了，比起看什么Legendre 16点逼近应该舒服一些。
                  </p>
                  <p>
                  </p>
                  <pre class="highlight "># 生成一个随机数
smm_one = function(m) max(abs(x &lt;- rnorm(m))) / sd(x)
# 生成n个随机数
rsmm = function(n, m) replicate(n, smm_one(m))
# 求分位数
qsmm = function(q, n = 10000, m) quantile(rsmm(n, m), q)
# 累积概率密度函数
psmm = function(x, n = 10000, m) mean(rsmm(n, m) &lt;= x)</pre>
                  <p>
                  </p>
                  <p>
                   话说你动作真够快的，呼哧一下就是一个R包上CRAN了，这种顺利程度建议出门买彩票去。
                  </p>
                 </div>
                 <!-- .bbp-reply-content -->
                </div>
                <!-- .reply -->
                <div class="bbp-reply-header" id="post-388129">
                 <div class="bbp-meta">
                  <span class="bbp-reply-post-date">
                   2013年12月22日 上午5:21
                  </span>
                  <span class="bbp-header">
                   回复：
                   <a class="bbp-topic-permalink" href="http://cos.name/cn/topic/141662/">
                    用knitr写报告百分号的意义冲突
                   </a>
                  </span>
                  <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/141662/#post-388129">
                   6 楼
                  </a>
                  <span class="bbp-admin-links">
                  </span>
                 </div>
                 <!-- .bbp-meta -->
                </div>
                <!-- #post-388129 -->
                <div class="even bbp-parent-forum-999 bbp-parent-topic-141662 bbp-reply-position-6 user-id-1 post-388129 reply type-reply status-publish hentry">
                 <div class="bbp-reply-author">
                  <a class="bbp-author-avatar" href="http://cos.name/cn/profile/1/" rel="nofollow" title="查看谢益辉的档案">
                   <img src="http://sdn.geekzu.org/avatar/1022d8e6ebc94e8f6bca9a86cebe312a?s=80&amp;d=monsterid&amp;r=g"/>
                  </a>
                  <br/>
                  <a class="bbp-author-name" href="http://cos.name/cn/profile/1/" rel="nofollow" title="查看谢益辉的档案">
                   谢益辉
                  </a>
                  <br/>
                  <div class="bbp-author-role">
                   站长
                  </div>
                 </div>
                 <!-- .bbp-reply-author -->
                 <div class="bbp-reply-content">
                  <p>
                   升级？
                  </p>
                 </div>
                 <!-- .bbp-reply-content -->
                </div>
                <!-- .reply -->
                <div class="bbp-reply-header" id="post-388128">
                 <div class="bbp-meta">
                  <span class="bbp-reply-post-date">
                   2013年12月22日 上午5:21
                  </span>
                  <span class="bbp-header">
                   回复：
                   <a class="bbp-topic-permalink" href="http://cos.name/cn/topic/144377/">
                    新发布了一个在线web scraping的站点,需要抓取网页数据的可以看看
                   </a>
                  </span>
                  <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/144377/#post-388128">
                   7 楼
                  </a>
                  <span class="bbp-admin-links">
                  </span>
                 </div>
                 <!-- .bbp-meta -->
                </div>
                <!-- #post-388128 -->
                <div class="odd bbp-parent-forum-999 bbp-parent-topic-144377 bbp-reply-position-7 user-id-1 post-388128 reply type-reply status-publish hentry">
                 <div class="bbp-reply-author">
                  <a class="bbp-author-avatar" href="http://cos.name/cn/profile/1/" rel="nofollow" title="查看谢益辉的档案">
                   <img src="http://sdn.geekzu.org/avatar/1022d8e6ebc94e8f6bca9a86cebe312a?s=80&amp;d=monsterid&amp;r=g"/>
                  </a>
                  <br/>
                  <a class="bbp-author-name" href="http://cos.name/cn/profile/1/" rel="nofollow" title="查看谢益辉的档案">
                   谢益辉
                  </a>
                  <br/>
                  <div class="bbp-author-role">
                   站长
                  </div>
                 </div>
                 <!-- .bbp-reply-author -->
                 <div class="bbp-reply-content">
                  <p>
                   这个太猛了。本来是要召唤土地公公，结果一顺手就招来神龙。
                  </p>
                 </div>
                 <!-- .bbp-reply-content -->
                </div>
                <!-- .reply -->
                <div class="bbp-reply-header" id="post-388124">
                 <div class="bbp-meta">
                  <span class="bbp-reply-post-date">
                   2013年12月22日 上午2:37
                  </span>
                  <span class="bbp-header">
                   回复：
                   <a class="bbp-topic-permalink" href="http://cos.name/cn/topic/140677/">
                    请问可以在knitr包的chunk set 区域中使用Sexpr{}吗？
                   </a>
                  </span>
                  <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/140677/#post-388124">
                   2 楼
                  </a>
                  <span class="bbp-admin-links">
                  </span>
                 </div>
                 <!-- .bbp-meta -->
                </div>
                <!-- #post-388124 -->
                <div class="even bbp-parent-forum-999 bbp-parent-topic-140677 bbp-reply-position-2 user-id-1 post-388124 reply type-reply status-publish hentry">
                 <div class="bbp-reply-author">
                  <a class="bbp-author-avatar" href="http://cos.name/cn/profile/1/" rel="nofollow" title="查看谢益辉的档案">
                   <img src="http://sdn.geekzu.org/avatar/1022d8e6ebc94e8f6bca9a86cebe312a?s=80&amp;d=monsterid&amp;r=g"/>
                  </a>
                  <br/>
                  <a class="bbp-author-name" href="http://cos.name/cn/profile/1/" rel="nofollow" title="查看谢益辉的档案">
                   谢益辉
                  </a>
                  <br/>
                  <div class="bbp-author-role">
                   站长
                  </div>
                 </div>
                 <!-- .bbp-reply-author -->
                 <div class="bbp-reply-content">
                  <p>
                  </p>
                  <pre class="highlight ">&lt;&lt;lsheng,fig.cap=paste('得分前', geshu, '的省')&gt;&gt;=
qplot(x=year,y=score,data=lsheng,color=diqu,geom='line')
@</pre>
                  <p>
                  </p>
                 </div>
                 <!-- .bbp-reply-content -->
                </div>
                <!-- .reply -->
                <div class="bbp-reply-header" id="post-388121">
                 <div class="bbp-meta">
                  <span class="bbp-reply-post-date">
                   2013年12月22日 上午2:17
                  </span>
                  <span class="bbp-header">
                   回复：
                   <a class="bbp-topic-permalink" href="http://cos.name/cn/topic/140952/">
                    吐槽一下R的正则表达式
                   </a>
                  </span>
                  <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/140952/#post-388121">
                   3 楼
                  </a>
                  <span class="bbp-admin-links">
                  </span>
                 </div>
                 <!-- .bbp-meta -->
                </div>
                <!-- #post-388121 -->
                <div class="odd bbp-parent-forum-999 bbp-parent-topic-140952 bbp-reply-position-3 user-id-1 post-388121 reply type-reply status-publish hentry">
                 <div class="bbp-reply-author">
                  <a class="bbp-author-avatar" href="http://cos.name/cn/profile/1/" rel="nofollow" title="查看谢益辉的档案">
                   <img src="http://sdn.geekzu.org/avatar/1022d8e6ebc94e8f6bca9a86cebe312a?s=80&amp;d=monsterid&amp;r=g"/>
                  </a>
                  <br/>
                  <a class="bbp-author-name" href="http://cos.name/cn/profile/1/" rel="nofollow" title="查看谢益辉的档案">
                   谢益辉
                  </a>
                  <br/>
                  <div class="bbp-author-role">
                   站长
                  </div>
                 </div>
                 <!-- .bbp-reply-author -->
                 <div class="bbp-reply-content">
                  <p>
                   用R的正则表达式时就是在三次元和二次元空间之间游离……有时候二次转义，有时候三次，有时候又不需要。有时候R core都能自己砸到自己的脚，比如这个存在了十几年的Sweave bug（中括号内大括号不需要转义）我吼了几嗓子之后Ripley大人明明修正了后来不知何故又滚回了：https://github.com/wch/r-source/commit/52b0a46e
                  </p>
                 </div>
                 <!-- .bbp-reply-content -->
                </div>
                <!-- .reply -->
                <div class="bbp-reply-header" id="post-388119">
                 <div class="bbp-meta">
                  <span class="bbp-reply-post-date">
                   2013年12月22日 上午2:03
                  </span>
                  <span class="bbp-header">
                   回复：
                   <a class="bbp-topic-permalink" href="http://cos.name/cn/topic/141392/">
                    R做回归怎么把小于2e-16的Pr(&gt;|t|)显示出来？
                   </a>
                  </span>
                  <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/141392/page/2/#post-388119">
                   20 楼
                  </a>
                  <span class="bbp-admin-links">
                  </span>
                 </div>
                 <!-- .bbp-meta -->
                </div>
                <!-- #post-388119 -->
                <div class="even bbp-parent-forum-999 bbp-parent-topic-141392 bbp-reply-position-20 user-id-1 post-388119 reply type-reply status-publish hentry">
                 <div class="bbp-reply-author">
                  <a class="bbp-author-avatar" href="http://cos.name/cn/profile/1/" rel="nofollow" title="查看谢益辉的档案">
                   <img src="http://sdn.geekzu.org/avatar/1022d8e6ebc94e8f6bca9a86cebe312a?s=80&amp;d=monsterid&amp;r=g"/>
                  </a>
                  <br/>
                  <a class="bbp-author-name" href="http://cos.name/cn/profile/1/" rel="nofollow" title="查看谢益辉的档案">
                   谢益辉
                  </a>
                  <br/>
                  <div class="bbp-author-role">
                   站长
                  </div>
                 </div>
                 <!-- .bbp-reply-author -->
                 <div class="bbp-reply-content">
                  <p>
                   如果楼主不是处女座的话，我就不能理解这个帖子的意思了。
                  </p>
                 </div>
                 <!-- .bbp-reply-content -->
                </div>
                <!-- .reply -->
                <div class="bbp-reply-header" id="post-388118">
                 <div class="bbp-meta">
                  <span class="bbp-reply-post-date">
                   2013年12月22日 上午2:02
                  </span>
                  <span class="bbp-header">
                   回复：
                   <a class="bbp-topic-permalink" href="http://cos.name/cn/topic/141401/">
                    请教一个时间序列数据的预测问题
                   </a>
                  </span>
                  <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/141401/#post-388118">
                   3 楼
                  </a>
                  <span class="bbp-admin-links">
                  </span>
                 </div>
                 <!-- .bbp-meta -->
                </div>
                <!-- #post-388118 -->
                <div class="odd bbp-parent-forum-1015 bbp-parent-topic-141401 bbp-reply-position-3 user-id-1 post-388118 reply type-reply status-publish hentry">
                 <div class="bbp-reply-author">
                  <a class="bbp-author-avatar" href="http://cos.name/cn/profile/1/" rel="nofollow" title="查看谢益辉的档案">
                   <img src="http://sdn.geekzu.org/avatar/1022d8e6ebc94e8f6bca9a86cebe312a?s=80&amp;d=monsterid&amp;r=g"/>
                  </a>
                  <br/>
                  <a class="bbp-author-name" href="http://cos.name/cn/profile/1/" rel="nofollow" title="查看谢益辉的档案">
                   谢益辉
                  </a>
                  <br/>
                  <div class="bbp-author-role">
                   站长
                  </div>
                 </div>
                 <!-- .bbp-reply-author -->
                 <div class="bbp-reply-content">
                  <p>
                   这个问题中日期真的有意义吗？我估计没有（当然你可以把日期加进来检验一下是否显著）。换句话说，只有某一天当中的观测时刻是有用的自变量，那么你把横坐标上的四条线压缩到同一天内做模型就好了（比如对时间回归之类的，不是所有时间序列非用ARIMA不可），预测的时候可以不必管具体日期。
                  </p>
                 </div>
                 <!-- .bbp-reply-content -->
                </div>
                <!-- .reply -->
                <div class="bbp-reply-header" id="post-388115">
                 <div class="bbp-meta">
                  <span class="bbp-reply-post-date">
                   2013年12月22日 上午1:50
                  </span>
                  <span class="bbp-header">
                   回复：
                   <a class="bbp-topic-permalink" href="http://cos.name/cn/topic/141637/">
                    回归分析中因变量需要服从正态分布嘛？
                   </a>
                  </span>
                  <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/141637/#post-388115">
                   8 楼
                  </a>
                  <span class="bbp-admin-links">
                  </span>
                 </div>
                 <!-- .bbp-meta -->
                </div>
                <!-- #post-388115 -->
                <div class="even bbp-parent-forum-991 bbp-parent-topic-141637 bbp-reply-position-8 user-id-1 post-388115 reply type-reply status-publish hentry">
                 <div class="bbp-reply-author">
                  <a class="bbp-author-avatar" href="http://cos.name/cn/profile/1/" rel="nofollow" title="查看谢益辉的档案">
                   <img src="http://sdn.geekzu.org/avatar/1022d8e6ebc94e8f6bca9a86cebe312a?s=80&amp;d=monsterid&amp;r=g"/>
                  </a>
                  <br/>
                  <a class="bbp-author-name" href="http://cos.name/cn/profile/1/" rel="nofollow" title="查看谢益辉的档案">
                   谢益辉
                  </a>
                  <br/>
                  <div class="bbp-author-role">
                   站长
                  </div>
                 </div>
                 <!-- .bbp-reply-author -->
                 <div class="bbp-reply-content">
                  <p>
                   回复 第6楼 的 nuomin：确切地说是Y|X服从正态分布，或者说误差项服从正态分布，而不是Y服从正态分布。Y本身完全可以不服从正态分布。
                  </p>
                 </div>
                 <!-- .bbp-reply-content -->
                </div>
                <!-- .reply -->
                <div class="bbp-reply-header" id="post-388114">
                 <div class="bbp-meta">
                  <span class="bbp-reply-post-date">
                   2013年12月22日 上午1:47
                  </span>
                  <span class="bbp-header">
                   回复：
                   <a class="bbp-topic-permalink" href="http://cos.name/cn/topic/141662/">
                    用knitr写报告百分号的意义冲突
                   </a>
                  </span>
                  <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/141662/#post-388114">
                   4 楼
                  </a>
                  <span class="bbp-admin-links">
                  </span>
                 </div>
                 <!-- .bbp-meta -->
                </div>
                <!-- #post-388114 -->
                <div class="odd bbp-parent-forum-999 bbp-parent-topic-141662 bbp-reply-position-4 user-id-1 post-388114 reply type-reply status-publish hentry">
                 <div class="bbp-reply-author">
                  <a class="bbp-author-avatar" href="http://cos.name/cn/profile/1/" rel="nofollow" title="查看谢益辉的档案">
                   <img src="http://sdn.geekzu.org/avatar/1022d8e6ebc94e8f6bca9a86cebe312a?s=80&amp;d=monsterid&amp;r=g"/>
                  </a>
                  <br/>
                  <a class="bbp-author-name" href="http://cos.name/cn/profile/1/" rel="nofollow" title="查看谢益辉的档案">
                   谢益辉
                  </a>
                  <br/>
                  <div class="bbp-author-role">
                   站长
                  </div>
                 </div>
                 <!-- .bbp-reply-author -->
                 <div class="bbp-reply-content">
                  <p>
                   <code>
                    "\\%"
                   </code>
                  </p>
                 </div>
                 <!-- .bbp-reply-content -->
                </div>
                <!-- .reply -->
               </li>
               <!-- .bbp-body -->
               <li class="bbp-footer">
                <div class="bbp-reply-author">
                 作者
                </div>
                <div class="bbp-reply-content">
                 帖子
                </div>
                <!-- .bbp-reply-content -->
               </li>
               <!-- .bbp-footer -->
              </ul>
              <!-- #topic-0-replies -->
              <div class="bbp-pagination">
               <div class="bbp-pagination-count">
                查看 15 个帖子 - 46 到 60（总计 8,682 个）
               </div>
               <div class="bbp-pagination-links">
                <a class="prev page-numbers" href="http://cos.name/cn/profile/1/replies/page/3/">
                 ←
                </a>
                <a class="page-numbers" href="http://cos.name/cn/profile/1/replies/">
                 1
                </a>
                <span class="page-numbers dots">
                 …
                </span>
                <a class="page-numbers" href="http://cos.name/cn/profile/1/replies/page/3/">
                 3
                </a>
                <span class="page-numbers current">
                 4
                </span>
                <a class="page-numbers" href="http://cos.name/cn/profile/1/replies/page/5/">
                 5
                </a>
                <span class="page-numbers dots">
                 …
                </span>
                <a class="page-numbers" href="http://cos.name/cn/profile/1/replies/page/579/">
                 579
                </a>
                <a class="next page-numbers" href="http://cos.name/cn/profile/1/replies/page/5/">
                 →
                </a>
               </div>
              </div>
             </div>
            </div>
            <!-- #bbp-user-replies-created -->
           </div>
          </div>
         </div>
        </div>
        <!-- .entry-content -->
        <footer class="entry-meta">
        </footer>
        <!-- .entry-meta -->
       </article>
       <!-- #post -->
      </div>
      <!-- #content -->
     </div>
     <!-- #primary -->
     <div class="sidebar-container" id="tertiary" role="complementary">
      <div class="sidebar-inner">
       <div class="widget-area">
        <aside class="widget bbp_widget_login" id="bbp_login_widget-2">
         <h3 class="widget-title">
          登录
         </h3>
         <form action="http://cos.name/cn/wp-login.php" class="bbp-login-form" method="post">
          <fieldset>
           <legend>
            登录
           </legend>
           <div class="bbp-username">
            <label for="user_login">
             用户名:
            </label>
           </div>
           <div class="bbp-password">
            <label for="user_pass">
             密码:
            </label>
           </div>
           <div class="bbp-remember-me">
            <label for="rememberme">
             记住用户名
            </label>
           </div>
           <div class="bbp-submit-wrapper">
            <button class="button submit user-submit" id="user-submit" name="user-submit" tabindex="104" type="submit">
             登录
            </button>
           </div>
           <div class="bbp-login-links">
            <a class="bbp-register-link" href="http://cos.name/cn/wp-login.php?action=register" title="注册">
             注册
            </a>
            <a class="bbp-lostpass-link" href="http://cos.name/cn/wp-login.php?action=lostpassword" title="忘记密码">
             忘记密码
            </a>
           </div>
          </fieldset>
         </form>
        </aside>
        <aside class="widget widget_text" id="text-7">
         <h3 class="widget-title">
          搜索
         </h3>
         <div class="textwidget">
          <form action="http://www.google.com/search" id="bbp-search-form" method="get" onsubmit="Gsitesearch(this)" role="search">
           <div>
           </div>
          </form>
          <form id="bbp-search-form-baidu" onsubmit="g(this)" role="search">
           <div>
           </div>
          </form>
         </div>
        </aside>
        <aside class="widget widget_text" id="text-2">
         <h3 class="widget-title">
          新鲜事
         </h3>
         <div class="textwidget">
          <ul>
           <li>
            <a href="http://cos.name/cn/topics/">
             最新帖子
            </a>
           </li>
           <li>
            <a href="http://cos.name/cn/view/popular/">
             最热门主题
            </a>
           </li>
           <li>
            <a href="http://cos.name/cn/view/no-replies/">
             消灭零回复
            </a>
           </li>
          </ul>
         </div>
        </aside>
        <aside class="widget widget_text" id="text-3">
         <h3 class="widget-title">
          RSS订阅
         </h3>
         <div class="textwidget">
          <ul>
           <li>
            <img src="http://cos.name/wp-includes/images/rss.png"/>
            <a href="http://cos.name/cn/topics/feed/">
             所有主题
            </a>
           </li>
           <li>
            <img src="http://cos.name/wp-includes/images/rss.png"/>
            <a href="http://cos.name/cn/forums/feed/">
             所有帖子
            </a>
           </li>
          </ul>
         </div>
        </aside>
       </div>
       <!-- .widget-area -->
      </div>
      <!-- .sidebar-inner -->
     </div>
     <!-- #tertiary -->
    </div>
    <!-- #main -->
    <footer class="site-footer" id="colophon" role="contentinfo">
     <div class="site-info">
      版权所有 © 2014 统计之都 | 由
      <a href="http://wordpress.org/">
       WordPress
      </a>
      构建 | 主题修改自
      <a href="http://wordpress.org/themes/twentythirteen">
       Twenty Thirteen
      </a>
     </div>
     <!-- .site-info -->
    </footer>
    <!-- #colophon -->
   </div>
   <!-- #page -->
  </div>
 </body>
</html>